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We consider the statistical mechanics of interfering transmissions in a wireless communications 
protocol. In this case, a connection between two nodes requires all other nodes within communication 
distance of the given two nodes to remain quiet on the given channel. This leads to an interesting 
problem of dimers on a lattice, with a restriction that no two dimers can overlap or be nearest 
neighbors. We consider both an equilibrium and a non-equilibrium, "greedy" dynamics for the links; 
the equilibrium properties of the model are found to exhibit an interesting spin glass transition at 
maximum density on certain lattices, while the greedy construction is related to the problem of 
random sequential adsorption. 



The study of networks has identified a variety of in- 
teresting networks with different topological properties 
[1], with recent interest considering dynamical processes 
such as flows on the networks [2]. Wireless communica- 
tion networks provide a practically important example 
of information flow on a network. In this paper, we con- 
sider the problem of the MAC (media access control) 
layer, which lies below the routing layer and allows dif- 
ferent nodes access to given frequencies. This layer must 
enable a large number of nodes to communicate over a fi- 
nite number of frequencies with finite signal range, while 
optimizing the latency, throughput, and fairness of the 
protocol [3]. In this paper, we show that even the sim- 
plest model of the MAC layer gives rise to rich behavior, 
including glassy and non-equilibrium dynamics. 

Recently, the problem of the MAC layer in such a pro- 
tocol has been cast as a matching problem [4] . Consider 
some set of nodes, each with given signal range. As a sim- 
ple approximation, the signal range leads to a directed 
network of possible communication links between nodes: 
there is a link from node i to node j if node j can receive 
a transmission from node %. As a further simplification, 
we consider the case in which the network is undirected, 
so that whenever node j can receive node i, then node i 
can also receive node j. 

Suppose there is only one allowed frequency. In this 
case, suppose two nodes, i and j establish a communica- 
tion link. Then, all nodes which are within range of node 
i or node j must remain silent, staying off the frequency 
in question, to avoid interference. Thus, it is possible for 
a given subset of links to be active simultaneously if and 
only if, for all nodes i, if node i is in an active link then 
i has exactly one neighbor j which is also active. Equiv- 
alently, this is the problem of finding a set of dimers in 
which nodes in different dimers are at least a distance 
of 2 from each other on the network; the presence of a 
dimer indicates an active link between the two nodes. 
The problem of multiple frequencies is similar, and leads 
to a coloring problem [5]. In Fig. 1, we show an exam- 
ple network and various possibilities of both allowed and 
disallowed sets of communication links. 



In this paper, we study the statistical mechanics of this 
problem. We consider two different dynamics to generate 
the links. In the first case, we consider summing over all 
possible allowed sets of dimers, weighting the covering 
by a fugacity for the number of dimers. By taking the 
fugacity to be large, we are thus able to find high density 
configurations of dimers on the network. We solve this 
problem when the network is a Bethe lattice. Interest- 
ingly, the results depend on the coordination number of 
the Bethe lattice: for coordination number greater than 
5, there appears to be a spin glass transition at a finite 
fugacity, while for coordination number less than 5, this 
transition is fixed at infinite fugacity. We also consider 
the problem on various finite dimensional lattices and 
speculate on the phase diagram. The problem is related 
to a recently considered statistical mechanics of dimers in 
the limit of infinite repulsion between neighboring dimers 
[6] , but we identify the novel possibility of the spin glass 
phase in this limit. 




FIG. 1. a) An example of an allowed configuration of 
dimers. Solid lines denote the presence of a dimer (active 
link), dashed lines denote a link on the network which has 
no dimer (an inactive link), b) A disallowed configuration of 
dimers. The two dimers are within a distance one of each 
other. 

The second dynamics we consider is a "greedy" dy- 
namics, in which links are added one at a time to the 
network, each link being added at random among the al- 
lowed possibilities. This dynamics is a version of random 
sequential adsorption [7]. We will see that the greedy 
dynamics does very poorly on graphs with high coordi- 
nation number. 

To motivate the equilibrium and greedy dynamics, we 
start with the following more general set of dynamics: 
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at every instant, wc suppose that some set of dimers is 
present. We then suppose that dimers can be created, 
when a link is activated, and disappear, when two nodes 
which were communicating stop using the channel. Both 
processes happen with given rates as follows. We suppose 
that at every instant of time each link which is inactive 
has some rate of activating (if that does not interfere with 
other links), and we set this rate to unity. We assign each 
dimer some rate of disappearing in any instant of time, 
and we set this rate to I /(f). Then, in the stationary dis- 
tribution of this dynamics, the probability of finding any 
given allowed configuration of dimers is Z^ 1 ^ , where 
Nd is the number of dimers in the given configuration 
and 



z = J2<p Nd (c), 



(i) 



where c is a configuration of dimers and Nd(c) is the 
number of dimers in the configuration. 

Then, the larger <f> is, the larger the average number 
of active links will be, thus increasing throughput. How- 
ever, a very large <f> means that the configuration of active 
links changes slowly in time; the system will tend to re- 
main stuck in a given configuration for a long time, which 
is also undesirable. We consider two cases below. First, 
the stationary dynamics at finite <f>. Next, we consider 
4> = oo, and start from an initial configuration with no 
links active; in this case, links can be activated but are 
never deactivated. This leads to a random sequential ad- 
sorption problem, and we will see that the final density 
of links in this case is far below the maximum possible in 
finite (j>. 

Equilibrium Dynamics on the Bethe Lattice — We start 
with the equilibrium case, solving for the partition func- 
tion (1). This problem is equivalent to that studied in [6] 
in the limit a — > — oo (in that paper's notation). How- 
ever, we consider the possibility of different spin glass 
and ordering transitions in this paper. 

We consider a Bethe lattice with fixed branching ratio, 
k, so that every node has q = k + 1 neighbors. There are 
a number of different possible boundary conditions con- 
sidered in the literature [8]. One possibility is to assume 
that the lattice indeed is a tree with given branching ratio 
and with some number of leaves; in this case, appropri- 
ate boundary conditions must be assigned at the leaves 
of the tree and one must consider only spins which are re- 
mote from the boundary. This procedure is reasonable in 
unfrustrated systems, and in general for systems outside 
of a spin glass phase. However, frustrated systems in the 
spin glass phase are sensitive to the boundary conditions 
which fix the degree of frustration [9]. Then, a better 
choice may be a random graph with fixed connectivity 
k + l. 

We will solve the problem with a set of recurrence equa- 
tions on the Bethe lattice, and identify the onset of a spin 
glass phase. In the spin glass phase, which we do not 



treat, the system becomes sensitive to boundary condi- 
tions. Outside the spin glass phase, the simple solution 
with recurrence equations is adequate to describe both 
the case of a tree with random boundary conditions and 
the case of a random graph. 

To obtain the recurrence equations, we consider a tree 
with one of three possible states for the root of the tree 
(the root has only k neighbors). First, the root can be 
connected to one of its daughters via a dimer (we call 
this active). Second, the root can be not connected to 
any of its daughters, but have at least one of its daugh- 
ters involved in an active link (we refer to this as locked) . 
Third, the root can be not connected to any of its daugh- 
ters, and have none of its daughters active either (we re- 
fer to this as free). These three possibilities are shown in 
Fig. 2(a-c). 
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FIG. 2. a) Active configuration: root connected to daugh- 
ter by an active link, b) Locked configuration: daughter in- 
volved in an active link, preventing root from forming a link 
to any other node, c) Free configuration: none of the daugh- 
ters have active links, d) Joining k + l trees as described in 
text. 

Let the partition function of the dimer problem on a 
tree in each of the three cases be Z\, Z2, Z 3 respectively. 
Then take k of these trees and join them together to 
make a tree of one higher level, and compute the resulting 
partition functions Z\,Z 2 , Z 3 for the larger tree with the 
three different boundary conditions for that tree. 

We have that Z\ = <pkZ 3 (Z 2 + Z 3 ) k ^ 1 , as one of the k 
daughters must be free to have a dimer with the new root, 
and all the other daughters must be locked or free (they 
cannot be active or the root could not make a connection 
to the free daughter). Then, Z 3 = (Z 2 + Z 3 ) k , as each of 
the k daughters of the new root must be either locked or 
free, but cannot be active, so each subtree has partition 
function Z 2 + Z 3 . Finally, Z 2 = {Z 1 + Z 2 + Z 3 ) k - (Z 2 + 
Z 3 ) k , as in this case at least one daughter must be active; 
so, we consider all possible states of the daughters, [Z\ + 
Z 2 + Z 3 ) k , and subtract off those states with no daughters 
active, Z 2 = (Z 2 + Z 3 ) k . 

It is convenient to rescale the partition functions Zi 
by a constant in looking for a fixed point. We choose to 
scale Zi — > Zi/Z 3 , so Z 3 = 1. The rescaled equations 
become: 
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Z 3 = l; Zi =#(Z 2 + l) fe -V(Z 2 + l) fc ; (2) 
Z 2 = (Zi + Z 2 + l) fc - (Z 2 + l) fc /(Z 2 + l) fc . 

We have numerically solved Eqs. (2) and found that for 
small 4> the solution converges to a fixed point, which we 
refer to as a liquid phase. For larger 0, there is a range 
for which the fixed point is stable, but there is also a sta- 
ble set of solutions which oscillate with period three. The 
oscillating solution corresponds to a solid phase with bro- 
ken translational symmetry, in which at one given level 
of the tree most of the nodes are free, at the next higher 
level most of the nodes are active to a daughter, and at 
the next higher level most of the nodes are locked. 

At even higher values of <f>, the stationary solution be- 
comes unstable and only the solid solution remains. The 
existence of the solid solutions depends on the boundary 
conditions on the tree. If the lattice actually is a balanced 
tree with given branching ratio and depth, then the sys- 
tem has a phase transition to a solid phase. However, on 
a system on a random graph, the presence of loops will 
frustrate the solid phase. Alternately, on a tree in which 
the boundary conditions on the leaves are chosen suffi- 
ciently randomly (for example, the tree is not balanced, 
but instead unbalanced, so that different leaves are at 
different depths below the root), again we can stabilize 
the liquid phase. For given 0, there is only one liquid 
phase (the possibility of a liquid-gas transition as in [6] 
does not happen in this system). 

At large </>, we can solve Eq. (2) in the liquid phase 
giving z[ iquid = fc(fc+i)/(2fc+i)^(fc+i)/(2fc+i)^<z«i<i _ 

k k/(2k+l) (j) k/(2k+l)^ Z li q uid = L To find probabn _ 

ity that a given site has a dimer leaving it, we must 
then connect q = k + 1 of these trees into 1 tree, as 
in Fig. 2(d), rather than k trees as above. After carrying 
this out, one finds that for large <f> a site has probability 
(k + l)/(2fc+ 1) - O^-VPfc+i)) of having an active link 
emanating from it. Defining N to be the total number 
of sites and p = 2Nd/N, we get a limiting dimer density 
at large <p of p = (k + l)/(2fc + 1) = q/(2q - 1). 

The limiting density in the liquid phase can be under- 
stood on physical grounds. For any lattice with fixed 
coordination number k + 1 and with no loops of length 3, 
the dimer density is at most (k + I) /(2k +1): consider a 
pair of sites 1, 2 connected by a dimer. Each of those two 
sites connects to k other sites, giving a total of 2k other 
sites connecting to the given two sites 1 and 2 (by assum- 
ing the absence of loops of length 3, the sites neighboring 
site 1 are distinct from those neighboring site 2). Thus, 
for each dimer, we can identify 2k sites which neighbor 
the given dimer; the total number of sites which neigh- 
bor dimers is then at least 2k /(k + 1): since each site has 
k + 1 neighbors, we may have overcounted the number 
of sites which neighbor dimers by a factor k + 1. Then, 
the total number of sites which have active links, plus 
the total number of sites which neighbor a site with an 
active link is at least N D [2 + 2k/(k + 1)] < N. Thus, 



p = 2N D /N < (k + l)/(2k + 1). To achieve this den- 
sity, it is necessary that if a node does not have an active 
link then all k + 1 of its neighbors have active links as in 
Fig. 3(a). While it appears that the system does achieve 
this maximum density as <j> — > oo from the results above, 
we now consider the possibility of a spin glass transition 
in the system which may prevent it from reaching the 
maximum density. 




FIG. 3. a) Configuration to obtain highest density on Bethe 
lattice, b) Highest possible density configuration on the 
square lattice. 



To identify the spin glass transition, we follow the pro- 
cedure used for the Ising spin glass [10,8]. Let each of 
the k trees that we join together in the procedure above 
have slight random variations in the values of Z\, Z 2 , Z 3 
about the liquid solution above, resulting from random- 
ness in the boundary conditions. If the liquid solution is 
stable against these random fluctuations, then the Bethe 
lattice and random graph will behave identically in the 
thermodynamic limit. Otherwise, we go to a spin glass 
phase, and the complete graph may have different be- 
havior. Fixing Z3 = 1 as above, we have two physical 
quantities, Z\ and Z 2 for each subtree. Suppose first 
that on each subtree each of these two values has the 
same slight fluctuation about the liquid solution above. 
We can compute the change in Z\ and Z 2 to linear or- 
der from Eq. (2). The resulting linear transformation is 
found to have eigenvalues 

A± = _k±^-4k) (3) 

at <f> — 00. Let these eigenvalues have left eigenvectors 
v ± . Following [10], we consider the mean-square fluctu- 
ations M ± = I E l ^ ± (^i ~ Z f quid )\ 1 . Here, the overline 
denotes averaging over the random fluctuations, while 
the sum over i = 1,2 projects (Zi — Z* qul ) onto the 
corresponding right eigenvector. If on each subtree the 
fluctuations about the mean are random, then the mean- 
square fluctuation at the next higher level is given by 
M ± = \X\ 2 /kM ± , where the reduction by k is caused by 
averaging k random numbers together. When |A| 2 = k, 
we have a spin glass transition. Now, for k < 4, the two 
eigenvalues are complex conjugates with absolute value 
\fk at 4> — 00, so that the spin glass transition occurs 
precisely at <f> = 00, with no transition at (f> < 00. 
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For k > 4, we have a spin glass transition at finite </>. 
The eigenvalues are both real in this case, with only one 
of them greater than \fk in absolute value; this eigenvalue 
is negative. The presence of only one such eigenvalue 
suggests that the phase transition to the spin glass phase 
might be in the same universality class as the Ising spin 
glass on the random graph. 

Equilibrium in Finite Dimensions — The treatment 
above is on the Bethe lattice. The properties of the sys- 
tem on finite dimensional systems are also of interest. 
On a square lattice, the system has a number of possible 
ordered structures. One can show that the pattern of 
Fig. 3(b) is the highest possible density. Thus, we con- 
jecture that as a function of <fi the system has a phase 
transition from a liquid state at low 4> to the ordered 
state of Fig. 3(b) at large <f>, with possible additional 
transitions to other ordered phases. A very interesting 
question is whether the system may have a glassy phase 
at some intermediate </>. 

Greedy Dynamics — The second problem we consider 
is the random sequential adsorption case at <j> = oo. Our 
technique follows that of [11] in a slightly different case 
of no interaction between dimers. We start at time t = 
with no links on the lattice. We then add allowed links at 
a unit rate per link. We consider the probability E m (t) 
that after time t a connected cluster of m sites has (1): no 
links on any of the to sites and (2): none of the neighbors 
of the given m sites has a link. These conditions imply 
that the to sites are free to participate in links as they are 
not blocked by any of their neighbors. It is interesting 
that the probability can be written as E m (t), indepen- 
dent of the structure of the connected cluster (there are 
many possible clusters with to sites). At t = 0, we have 
E m (t) = 1 for all to. At later times, one can derive the 
following differential equation for to > 0: 

d t E m - -(m - l)E m - [{k - l)m + 2}E m+1 - (4) 

[(k - l)m + 2}kE m+2 . 

There first term corresponds to the possibility of adding 
a dimer to any of the (to — 1) bonds that connect two 
sites, both in the cluster. The second term corresponds 
to the possibility of connecting a dimer to any of the 
[(k — l)m + 2] bonds which connect a site in the cluster to 
a site outside the cluster, while the last term corresponds 
to the possibility that a site which neighbors the cluster 
will connect to one of the k sites outside the cluster which 
neighbor the given site. 

The ansatz E m (t) = c^o-^)™" 1 solves Eq. (4) with 
d t a = -a - (k - l)a 2 - (k - l)fccr 3 and d t c = -(k + 
1)<jc — (k + l)k<7 2 c and initial conditions c(0) = <t(0) = 1. 

Define p(t) = 2Nu(t)/N to be the fraction of sites with 
a dimer emanating from them. We find dtp = (k + 1)E2- 

Numerical solution of the differential equation gives 
p(t — > oo) s=a .313 for k = 4. This is noticeably worse than 
the limiting value of p = 5/9 computed in the equilibrium 



case above. At larger k, the non-equilibrium dynamics 
does even worse. The differential equations can be solved 
asymptotically at large k, where for c, a » 1/k the 
equations simplify to dt<J = —k 2 o- 3 ,dtc — —k 2 a 2 c. The 
result for large k is 

p(t->oo)~ln(fc)/fc. (5) 

Discussion — We have studied a simple model inspired 
by recent developments in wireless communication. This 
model turns out to be a model of interacting dimers 
very similar to other models studied in statistical physics. 
From the point of view of statistical physics, one of the 
most interesting features of the model is the spin glass 
transition at <f> = oo for k < 4. In the language of statis- 
tical physics, this implies that the transition happens at 
zero temperature in the mean-field limit. This may sim- 
plify the treatment of this transition and may be useful 
for studying finite dimensional spin glasses. 

The spin glass transition is of interest for practical pur- 
poses also. In in the liquid phase heuristic algorithms 
will likely do a good job of rapidly finding allowed dimer 
configurations of given density, but in the glass phase it 
will be much more difficult to find good configurations 
and the system will tend to get stuck in certain configu- 
rations. Finally, the behavior of this system on regular 
lattices and random geometric graphs [12] in finite di- 
mensions may be interesting. 
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